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ABSTRACT 


An  approximate  method  has  been  devised  for  determining  the 
nonstationary  air-loads  on  an  elastic  wing  with  supersonic  edges  upon 
which  is  mounted  an  axially  symmetric  body.  The  flow  field  due  to  an 
oscillating  axially  symmetric  body  isolated  in  a supersonic  stream  is 
superimposed  on  the  flow  field  due  to  an  oscillating  elastic  wing  with 
supersonic  edges.  An  additional  velocity  potential  satisfying  the  basic 
linearized  differential  equation  of  fluid  motion  is  constructed  so  that 
the  required  boundary  conditions  are  satisfied  on  the  wing.  The  result 
is  a modification  of  the  original  source  strength  distribution  on  the 
isolated  wing.  This  modified  source  strength  distribution  is  calculated 
by  means  of  a numerical  integration  procedure.  The  nonstationary  air-loads 
for  the  wing  may  then  be  calculated  as  in  Reference  1. 
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INTRQDDCTICR 


Tne  trend  in  recent  high  speed  airplane  and  missile  designs 
has  been  toward  low  aspect  ratio  lifting  surfaces  mounted  on  bodies 
whose  diameters  are  large  fractions  of  the  total  lifting  surface 
span.  In  the  case  of  higher  aspect  ratio  configurations,  it  is  common 
practice  for  flutter  analysts  to  ignore  the  presence  of  the  body  or  to 
assume  that  the  wing  continues  uninterrupted  to  the  center  line  of  the 
body  when  calculating  unsteady  air-load  distributions  on  the  wing. 

For  moderate  aspect  ratios,  flutter  analysts  sometimes  treat  the  body 
as  an  infinite  reflecting  plane.  As  the  aspect  ratio  is  still  further 
reduced  so  that  the  body  diameter  becomes  proportionately  larger,  one 
is  led  to  ask  whether  the  effect  of  wing-body  interference  can  still 
be  ignored  in  computing  non-stationary  air-loads.  Certainly,  sufficient 
evidence  is  not  available  to  show  that  such  effects  are  unimportant.  It 
has  already  been  shown  that  the  presence  of  a body  near  a low  aspect 
ratio  lifting  surface  in  a steady  compressible  flow  appreciably  affects 
the  load  distribution  of  the  lifting  surface.  (See,  for  example,  Reference 
2).  One  is  therefore  inclined  to  suspect  that  wing-body  .interference  in 
non-stationary  flow  is  likewise  not  negligible.  Because  large  body 
diameter- wing  span  ratio  configurations  are  most  often  encountered  in 
high  speed  aircraft  and  missile  designs,  it  seems  appropriate  to  confine 
consideration  to  the  case  of  supersonic  flow. 

Since  the  linearized  problem  of  an  isolated  wing  oscillating 
in  a supersonic  stream  has  not  yet  been  solved  exactly,  it  cannot  be 
expected  that  an  exact  solution  for  the  ncn -stationary  supersonic  wing- 
body  problem  can  be  obtained  at  this  time.  However,  in  what  follows, 
an  approximate  method  will  be  derived  for  determining  the  effect  of  an 
axially  symmetric  fuselage  on  the  air-loads  acting  on  an  elastic  wing 
with  supersonic  edges  mounted  on  the  body  in  the  plane  of  the  body 
diameter. 


The  problem  of  the  body  alone  oscillating  in  a supersonic 
stream  is  treated  in  Reference  3.  Recently,  a group  at  Republic  Aviation 
Corporation  has  developed  an  approximate  method  for  computing  oscillating 
air-loads  bn  an  isolated  wing,  all  of  whose  edges  are  supersonic  (Reference 
1).  The  method  can  be  applied  to  an  elastic  wing  undergoing  arbitrary, 
harmonic  distortions.  The  formulation  of  the  present  problem  is  based 
on  the  assumptions  that  (a)  the  perturbations  from  a uniform  flow  are 
sufficiently  small  so  that  the  flow  equations  may  be  linearized,  and  (b) 
the  flow  field  around  the  wing  is  adequately  described  by  a linear  super- 
position of  the  velocity  potential  due  to  the  wing  alone,  the  velocity 
potential  due  to  the  body  alone,  and  an  additional  velocity  potential, 
which  when  added  to  the  first  two,  causes  the  boundary  conditions  on  the 
wing  to  be  satisfied.  The  velocity  potential  due  to  the  wing  alone 
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will  be  taken  from  the  work  of  Garrick  (Reference  1*),  the  potential 
for  the  body  alone  from  (Reference  5)$  the  wing-body  interference 
potential  is  developed  herein.  The  method  of  Reference  1 is  used  to 
obtain  numerical  solutions  of  the  equations. 
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SECTION  I 


CONSTRUCTION  OF  THE  SOLUTION 


The  velocity  potential  for  the  nonstationary  flow  over  a 
wing  with  supersonic  edges  has  been  derived  (under  the  assumption  of 
linearized  theory)  in  Reference  4.  In  that  derivation,  it  was  assumed 
that  the  actual  wing  may  be  replaced  by  a distribution  of  sources  and 
sinks  in  the  plane  z - 0.  It  was  shown  that  the  strength  of  a source 
at  any  point  on  the  wing  is  proportional  to  the  downwash  ( w*  ) at 
that  point.  Thus  for  harmonic  oscillations,  the  velocity  potential 
at  a point  (x,  y)  of  the  wing  for  the  flow  around  the  wing  alone  may 
be  written  as  . _ . 

[ur’r(Me  c°s(«Q  JfJf 

0 R 

where  u>  is  the  frequency  parameter  defined  as 

i3  the  frequency  of  oscillation,  M is  the  free  stream  Mach  number, 
c is  the  free  stream  speed  of  sound,  & is  the  Prandtl-Glauert  factor 


(1) 


(2) 


p 


(3) 


R may  be  written  as 


ft** 


U) 


and  D is  the  area  within  the  fore  Mach  cone  emanating  from  the  point 

(*»  y). 


The  velocity  potential  for  the  nonstationaiy  flow  about  a 
slender  body  of  revolution  has  been  derived  (under  the  assumptions  of 
linearized  theory)  in  Reference  3*  In  this  derivation  it  was  assumed 
that  the  bocfcr  may  be  replaced  by  a distribution  of  doublets  along  the 
axis  of  the  bocfy.  The  doublet  intensity  distribution  was  determined  by 
satisfying  the  condition  of  no  nonnal  velocity  at  the  surface  of  the 
body.  It  was  thereby  found  that  the  doublet  intensity  at  any  point 
is  proportional  to  the  product  of  the  body  cross-sectional  area(S)  at  that 
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point  times  the  vertical  velocity  of  the  body  axis  ( ) at  the  same 

point.  Thus  for  harmonic  oscillations,  the  velocity  potential  at  any 
point  in  space  for  the  flow  around  the  body  alone  may  be  written  (in 
cylindrical  coordinates,  (see  Figure  1)  as 


.if  (?) s '(%) f G) S(!)](*- Ve  cos(%P) 

r)  a 


(5) 


.igj  w.fflS<TX4r-»e'“*  cosffiP)  ^ J 

where  a prime  denotes  differentiation  with  respect  to  ^ and 

^skx-^)z-plrli 


(6) 


It  is  now  assumed  that  the  total  velocity  potential  •f  for  the 
wing-body  combination  may  be  expressed  as  the  sum  of  these  two  potentials 
plus  a correction  potential  % which  will  cause  the  boundary  condition 
to  be  satisfied  on  the  surface  of  the  wingj  that  is: 


(7) 


The  boundary  condition  on  the  wing  is 


5 If 

2Z 


_ 3* 


ur 


3 i 


'l*o 


Mi 

32 
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+ Mi 

32 


= ur 


AT 


2*0 


= oir 


ur 


(S) 


•2=0 
However, 

I 

32  U--0 

Since  the  velocity  potential  for  the  wing  alone  already  satisfies  the 
boundary  condition  for  the  motion  of  the  wing,  it  is  only  necessary  that 


(9) 
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*4  be  of  a form  such  that  the  vertical  velocity  that  it  induces 
at  a given  point.  on  the  wing  just  eancsls  the  downwash  induced  by  the 
body  at  the  same  point}  that  is, 


'2<Pd\  __ 


2-o 


(10) 


i-o 


An  alternate  formulation  of  this  problem  is  to  assume  that 
the  down-id  ah  produced  by  the  sources  representing  the  wing  must  be  modi- 
fied to  account  for  the  downwash  induced  by  the  body  so  that  the  boundary 
condition  on  the  wing  is  satisfied.  An  examination  of  equation  (5) 
shows  that  f b Is  always  zero  in  the  plane  z • 0 (0  E).  There- 

fore, the  potential  for  the  flow  over  the  wing  in  combination  with  the 
bocty-  is  of  the  same  form  as  the  potential  for  the  wing  alone,  but  the 
source  strength  distribution  is  modified  to  account  for  the  downwash  due 
to  the  bo<fy.  Thus, 

■■  t j fi. 


This  formulation  will  be  used  in  the  subsequent  development. 

Since  is  expressed  in  cylindrical  coordinates  and  it  is 
required  to  evaluate  at  z » 0.  is  written  as  follows : 

3*  3Z 


_ 3t- 

+ 39 

32 

, *30  di 

1-0 

(12) 


Noting  that 

.A  t x. 

r - y 


1*o 


(13) 


SL“  *1:° 

Since  remains  finite  on  the  wing,  the  first  term  of  equation  (12) 

vanishes.  Similarly, 


B-  tan  ' £ 


(15) 


1)  This  assumption  restricts  consideration  to  only  symmetrical  wing- 
body  configuration}  i.e.,  where  the  axis  of  the  body  lies  in  the  plane 
of  the  wing. 
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(16) 


and 


3Q  jr  _ L 

3*  i. o'  * 

Inspection  of  equation  (5)  shows  that  a f vnction’Vtf  (x, 
may  be  introducer!  such  that 


Vk cos  9[Vk(*,rjtj] 

and  therefore 

Ig  — suid  [*(*,#)] 

3 0 

At  z - 0,  Qn  and  r - yj  hence, 

!&— v*f**<; 


and 


3^j 

32 


where 


rL= 


iutt 


•f*o 
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‘4('sis('0e*tu  <jr"% /« 

- x d ? 

o P 

-&J  ^(t)$(&(x-%)<ii'orX’Vc os(%P) 


and  P is  now 


P=  JcZWfF 
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r , t ) 

(17) 

(18) 

(19) 

(20) 

(21) 

(22) 


i 


Substit,  \ting  equation  (20),  (21)  anrf  (22)  into  equation  (11) 
gives  the  complet  xpression  for  the  potential  at  point  (x,y)  on  the 
Td.ng  of  the  wing-b. _y  combination.  The  pressure  jump  at  this  point  i3 
then  obtained  from  the  potential  as 


The  lift  (positive  doim)  on  a streamwise  strip  of  the  wing  is  given  by 


and  the  stalling  moment  about  a spamrise  axis,  x0  , is 

Ms  =-  f 

S 


(22*) 


(25) 
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SECTION  II 


METHOD  OF  APPLICATION 


Up  to  the  present  time,  completely  analytical  evaluations 
of  integrals  of  the  form  of  equations  (1)  and  (21)  have  not  been 
accomplished.  For  rigid  bocfy  motions,  equation  (1)  has  been  approxi- 
mately evaluated  by  expanding  the  int'^rand  in  powers  of  uj  as  in 
Reference  5*  Ch  the  other  hand,  equation  (1)  has  been  approximately 
evalvated  even  for  the  case  of  arbitrary  wing  elastic  deformation 
modes  in  Reference  1 by  means  of  a numerical  integration  technique. 

In  this  paper,  both  of  these  methods  will  be  combined  to  evaluate 
equation  (21).  The  results  of  this  evaluation  may  then  be  used  directly 
in  calculating  the  nonsteady  air-loads  acting  on  the  wing  of  the 
wing-body  combination. 

The  first  step  in  this  technique  is  to  subdivide  the  wing 
into  a number  of  small  rectangular  boxes  of  dimensions  hx  and  h , 
as  shown  in  Figure  1.  Boxes  along  the  wing  trailing  edge  am  ma&e 
half -size  in  the  stream  direction,  in  accordance  with  the  procedure 
of  Reference  1.  The  body  is  also  divided  into  sections,  the  chordwise 
dimensions  of  which  are  made  equal  to  the  chordwise  dimension  of  the 
wing  boxes,  . The  size  of  the  boxes  is  not  important  from  the 

theoretical  point  of  view,  as  long  as  they  do  not  become  too  large 
in  the  stream  direction.  However,  certain  ratios  of  length  to  width 
of  the  rectangle  are  more  advantageous  for  computing  purposes.  (For 
a more  detailed  discussion  of  box  size  see  Reference  1). 

The  essence  of  the  method  of  Reference  1 is  that,  by  making  the 
boxes  small  enough,  the  downwash  over  any  box  can  be  considered  as 
essentially  constant  and  it  can  be  removed  from  within  the  integral 
of  equation  (1).  Then  by  expanding  the  numerator  of  the  remaining 
integrand  in  a power  series  in  co  , retaining  as  many  terms  as  are 
required  for  accuracy  in  practical  problems,  and  using  a mean  value 
theorem  where  the  power  series  expansion  is  not  needed,  the  double 
integral  can  be  evaluated  for  all  boxes  within  the  region  of  integration. 
This  procedure  will  be  adopted  for  the  evaluation  of  equation  (21). 

EVALUATION  OF  THE  BODY-INDUCED  DOWNWASH 

If,  then,  the  grid  of  boxes  is  made  fine  enough,  it  is 
reasonable  to  assume  that  w^  and  S are  constant  along  each  segment 
into  which  the  body  axis  has  been  divided.  Then,  the  values  of 
w^  and  S and  their  derivatives  with  respect  to  £j  at  any  point 
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along  a segment  are  taken  to  be  the  values  at  the  center  of  the 
segment*  Equation  (21)  then  becomes 

fc 

^ \ -iZo(x-f) 

Sj+iuunjsm  cx-f)e  c«s 
* 'J  p 

■where  the  symbol  ^ stands  for  summation  over  those  segments  of  the 
fuselage  centerline  included  within  the  fore  Mach  cone  emanating 
f ran  the  point  (x,  7)  on  the  wing,  and  the  symbol  stands  for 

integration  in  the  x direction  over  the  length  of  each  of  these 
segments.  For  a body  section  whose  center  line  is  cut  by  the  fore 
Mach  oone,  the  integration  is  performed  only  up  to  the  Mach  line. 
Thus  equation  (26)  may  be  written  as  follows: 


&*,?)*- 


( 


(27) 


For  the  case  of  segments  wherein  the  body  axis  is  not  cut  by  the  fore 
Mach  cone, 

e sin($P)J  f (2g) 


v(  < 

J *>, 


and 


i (*-f)  fULp) 

q-f)e  7 cos  ( Af 


Vjf/ 

* f-i 

Xj_|  ■ h^  « For  the  case  of  the  segnent  wherein  the  body 
is  cut  by  tne  fore  Mach  line, 

1,i  J e (30) 


where 

axis 
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and 


i 

/• 

*j-l 


-n  % -c-5 ->(*-$) 

Xx-0e 


cos 


(31) 


The  integrals  defined  by  equations  (28),  (29),  (30),  and 
(31)  may  be  evaluated  as  follows.  For  those  body  segments  which  are 
within  one  and  a half  box  lengths  in  the  flow  direction  from  the 
point  (x,y), (i.e. ,(x-  Xj_,)  ^ 3 </z  hx  ) » the  integrands  in  these 

integrals  may  be  expanded  in  powers  of  to  . Formulas  are  presented 
below  for  the  values  of  these  integrals  wherein  the  integrands  were 
expanded  to  the  second  power  of  to  . These  formulas  should  provide 
engineering  accuracy  for  most  problems.  Should  increased  accuracy  be 
desired,  it  can  be  achieved  by  using  these  formulas  with  a very  fine 
grid,  or  a coarse  _grid  and  more  accurate  formulas  which  include 
higher  powers  of  to  . Thus  for  a "near"  body  segment  not  cut  by  a 
Mach  line, 

I t - A,u>  + Az  u>* 

Zz*  60+£>,&  + ^ 4>z 


(32) 

(33) 


where 

al-  £ ft (55) 


6o-  la -*;-,AjSy  - J(r-Kj) 


Kff 


(36) 
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~ z fcx~  + & %flxcosh~l  (*t*0 

I ft 


(37) 


- (* - Xj-,)  Xj-,)2-^^ l-  ' ££i *td 

*ia&{0  +M  vm-j  VV7*-  (""l)[(*-  *h  t-f  yfi 

fa-^Wl}  C5S) 

Formulas  for  1^  and  Ig  for  the  case  of  "near"  segments 

wherein  the  body  axis  is  cut  by  the  Mach  line  can  be  obtained  from  th9 
above  formulas  by  putting  Xj  * x - ^ 

For  body  segments  which  are  more  than  one  and  a half  box 
lengths  removed  from  the  wing  point  (x,y),(x  - x . ^ 3/2  hx)  , the  ex- 
ponential and  trigonometric  terms  of  the  integrands  of  equations  (28), 
(29),  (30)  and  (31)  vary  very  slowly  over  the  length  of  the  segment 
provided  h_  is  reasonably  small.  Hence,  good  approximations  to  the 
values  of  these  integrals  may  be  obtained  by  considering  these  terms 
as  constant  and  equal  to  their  values  at  the  mid-point  of  the  segner.t 
over  which  the  integration  is  being  performed,  and  then  integrating 
the  remaining  part  of  the  integrand  directly.  That  is,  it  is  assumed 
for  "far"  segments  that 


r „ . 

s/'t, 


and 


e cos 


PCfJ 


(39) 


(ho) 


where 

v — Hi  £ SJj 
y*  “ z 


(hi) 
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For  those  body  segnents  wherein  the  body  center  line  is  not  cut  by  a 
Mach  line,  these  integrals  can  be  shorn  to  have  the  values 


(42) 


and 


(45) 


For  those  body  segments  wherein  the  body  center  line  is  cut  by 
a Mach  line,  the  same  equations  can  be  used  with  x^  ■ x -f»  * 

These  equations  make  it  possible  to  calculate  the  increments 
in  downwash  on  the  wing  due  to  the  body.  It  should  be  noted  that  these 
incremental  downwash  velocities  have  components  both  in  phase  and  out 
of  phase  with  the  wing  motion. 


DETERMINATION  OF  THE  NCN  STATIONARY  LIFT 
AND  MOMENT  DISTRIBUTIONS 


The  application  of  the  method  described  in  the  previous 
sections  to  the  calculation  of  the  unsteady  air-loads  on  the  wing 
of  a wing-body  combination  will  now  be  described.  Although  the  method 
is  generally  applicable  to  an  elastic  wing  experiencing  harmonic  os- 
cillations of  arbitrary  mode  shape,  for  simplicity  the  following 
derivation  will  be  restricted  to  a wing-body  configuration  undergoing 
rigid-body  symmetrical  motions.  In  order  to  determine  the  potential 
at  a point  (x,y)  on  the  wing,  equation  (11)  must  be  evaluated  over 
the  region  of  integration  D bounded  by  the  forward  Mach  cone  emanating 
from  (x,y),  the  wing  leading  edges,  and  the  sides  of  the  body.  The 
boxes  on  the  wing  which  contribute  to  ( x,y,t)  are  indicated  by  the 
cross-hatched  area  in  Figure  2.  In  keeping  with  the  assumptions  made 
previously,  the  downwash  over  a box  will  be  considered  constant  and 
equal  to  its  value  at  the  center  of  the  box.  Thus  the  geanetrical 
davnwash  due  to  the  motion  of  the  wing  is 


ur*r(%  %)  * * Vp 
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where  q^  and  are  the  generalized  coordinates  in  the  vertical 

translation  and  pitching  modes  respectively;  (4^  are  the 

coordinates  of  the  cenuroid  of  the  box  on  the  wing,  and  x0  is  the 
location  of  the  pitch  axis* 

The  downwash  at  ( df t » % ) induced  by  the  body  must  now 
be  determined.  A Mach  foreoone  is  drawn  from  (Sf- );  those  body 
segments  included  within  this  cone  (as  shown  in  Figure  2)  will 
contribute  to  the  induced  downwash  at  this  point. 

The  geometrical  downwash  at  any  point  along  the  body  center 
line  can  be  written  in  terms  of  the  generalized  coordinates  as  follows: 

* “off,  + * Vft 

From  this  one  can  immediately  write: 

* “n  (r*j)  = * <»[$<  us) 

Urkj'a  iWJ^A  (46) 


where  is  the  coordinate  of  the  center  of  the  jth  body  segment 
and  is  defined  by  equation  (i+1). 


Using  the  above  expressions  for  the  geometrical  downwash  of 
the  body,  one  can  rewrite  equations  (20)  and  (27)  as  follows  (dropping 
the  el“*  terms  for  convenience) 


where 

$ 


m 

(49) 


(50) 
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(51) 


Thus,  the  total  downwash  which  enters  into  equation  (11)  is 


(52) 


Equation  (11)  can  be  evaluated  in  accordance  with  the  "box" 
method  of  Reference  1,  leading  to  the  following  general  form  for  the 
velocity  potential  at  point  (x,y),  which  is  now  assumed  to  be  at  the 
center  of  box  "i": 


*(*,$)  - <fi 

= £*'"£  fy  hi '] 1<  *[iu>  ? + 


where 

• the  potential  at  the  center  of  box  i due  to  a constant  unit 
dcwnwash  velocity  at  box  j, 


y denotes  summation  over  all  boxes  included  within  the  region  of 
integration  D,  defined  previously.  Using  equations  (23)  and  (24), 
together  with  the  assumption  that  the  potential  at  the  center  of  a 
box  can  be  considered  as  the  mean  value  over  the  entire  box.  Reference 
1 shows  that  the  downward  load  acting  on  a chordwise  strip  of  the  wing 
is  given  by 

Vi  A*i+Zf>MrC  (57) 
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■where  the  index  i ranges  over  all  boxes  in  the  strip  from  leading  edge 
(LE)  to  trailing  edge  (TE)  and  is  the  value  of  the  ootential  at 

TE 

tne  trailing  edge  of  the  strip.  A X;  is  the  dimension  of  box  i in  the 
stream  direction.  Substituting  equation  (53)  into  equation  (57)  and 
non-dimensionalizing  yields  the  following  expressions  for  the  non- 
stationary sectional  lift  coefficients: 


where 


(58) 


L‘ = 4 jfcj  ?Vr,y  A/' 

-zT»  C 

bus 

with  k ■ — — — , the  local  reduced  frequency,  at  the  wing  root. 


(59) 

(60) 


In  a similar  manner,  using  equations  (25),  a general  expression 
for  the  pitching  moment  on  a strip  about  the  axis  of  pitch  xe  and 
equations  for  the  moment  coefficients  are  obtained. 


h 6Xi  *2^r£t/>T£  - z, y$  f.  A X: 
(M>  * L f (M3  + iM+)fx 

M ♦ ‘ M,  = ~ y (l,  * i Lt)  +J-p  p Axt$ 


(oi) 

(62) 


(63) 
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y_.  • ft;  / * ft; 

• , Uj  y , ft) 

'riaf  *“*f  v*  - JS*F«v‘< 


ft; 
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SECTION  III 


APPLICATION  TO  A SPECIFIC  PROBLEM 


The  application  of  the  method  described  in  the  previous 
sections  to  the  calculation  of  the  unsteady  air-loads  on  the  wing 
of  a wing-body  combination  will  now  be  illustrated.  Since  the 
purpose  of  this  example  will  he  merely  to  illustrate  the  application 
of  the  method  and  to  determine  the  order  of  magnitude  of  the 
effect  of  the  body  on  the  wing  air-loads,  a very  coarse  grid  will 
be  used  in  order  to  shorten  the  computations.  Rigid  body  translation 
and  pitching  motions  will  be  assumed  so  that  the  results  may 
be  compared  with  previously  obtained  results  for  the  wing  alone. 
Calculations  will  be  done  for  both  the  wing  alone  and  for  the 
wing-body  combination  so  that  the  effect  of  the  body  can  be 
determined  directly.  The  wing-alone  configuration  will  consist  of 
one  wing  of  the  wing-body  combination  reflected  about  the  wing 
root  chord  line.  The  wing  planform  chosen  is  the  same  as  the  one 
that  was  used  in  References  1 and  5»  This  will  also  allow  a com- 
parison to  be  made  between  the  results  obtained  with  a coarse  grid 
boxes;  Figure  2)  and  the  fine  grid  used  in  Reference  1 for  the 
wing  alone  case  (28  boxes;  Figure  1). 


Since  the  calculations  will  be  done  for  only  one  Mach 
number,  it  will  be  convenient  to  let  the  sides  of  the  rectangles  be 
in  the  ratio. 


<65) 


where  h and  h are  box  lengths  in  the  x and  y directions  respec- 
x y 

tively.  The  geometry  of  the  network,  therefore,  is  such  that  the 
body  axis  is  cut  by  a Mach  line  only  at  the  boundary  of  a body 
segnent. 


b 


The  following  parameters  are  used  in  this  analysis: 


M2  - 1.75, 


co  b 
V 


k **  .Oil 


iO 

co  M 
cp* 

xo 

- 2.25 

a 

- 1.75 

2bcoM2 


.18667 
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25 


Maximum  body  diameter 
Maximum  wing  span 

o 

Wing  leading  edge  sweep  angle  * j>0 

The  equation  for  the  radius  of  the  body  is  taken  to  be 
y - .U9l+S7x  - .053S63X5  (66) 

In  this  example  the  lift  and  moment  of  the  wing  are 
calculated  for  unit  values  of  — and  q . 

A detailed  drawing  of  the  wing  and  body  with  the  grid  of 
boxes  superimposed  is  shown  in  Figure  2. 

The  value  of  the  potential  for  boxes  close  to  the  wing 
leading  edge  is  obtained  by  weighting  the  result  obtained  for  the 
whole  box  in  accordance  with  the  percentage  of  area  of  the  whole 
box  that  is  contained  between  the  fore  Mach  cone  and  the  leading 
edge. 

1) 

The  results  of  calculations  using  equations  (57) 
through  (62+)  for  the  example  described  above  are  sumnarized  in 
Figures  3 through  12. 


T)  Alter  the  calculations  described  herein  were  completed,  the 
authors  of  Reference  1 devised  a scheme  whereby  the  lift  and 
moment  could  be  calculated  somewhat  more  simply  using  pressure 
coefficients  rather  than  velocity  potentials.  Since  the  principles 
involved  are  the  same,  the  calculations  are  presented  here  in  terms 
of  potentials. 
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SECTION  IV 


DISCUSSION  OF  RESULTS 


Formulas  have  been  derived  to  determine  the  effect  of  the 
fuselage  of  a wir^-body  combination  on  the  air-loads  anting  on  a 
■wing  in  nonstationary  purely  supersonic  flow*  The  equations  do 
not  readily  lend  themselves  to  analytic  evaluation.  However,  a 
procedure  has  been  presented  which  allows  a numerical  solution  of 
the  problem. 

Examination  of  equation  (5)  or  (21)  leads  to  the  conclusion 
that  the  most  accurate  rcsuDts  are  obtained  when  the  ratio  of  body 
diameter  to  wing  span  is  large  (providing  that  the  body  fineness 
ratio  is  large  so  that  linearized  theory  applies).  This  is  deduced 
from  the  fact  that  the  downwash  induced  by  the  body  has  a 1/r  type 
singularity  at  the  body  axis.  When  the  body  diameter  is  large,  no 
point  on  the  wing  is  very  close  to  the  body  axis  and  hence  the 
mathematical  singularity  which  does  not  actually  exist  physically 
need  not  be  dealt  with.  Practically,  the  presence  of  a singularity 
at  the  body  axis  is  not  really  a deterrent  to  the  use  of  this 
method  since  corrections  for  the  presence  of  the  body  are  probably 
not  required  when  the  tody  diameter  is  small  compared  to  the  wing 
span,  and  also  because  accuracy  is  more  important  in  the  region 
of  the  wing  tips  than  near  the  body. 

The  wing  lift  and  moment  due  to  harmonic  rigid  body  trans- 
lation and  pitching  about  the  mid-point  of  the  root  chord  for  the 
configuration  shown  in  Figure  2 are  presented  in  Figures  3 through 
10.  Curves  of  the  variation  of  the  unsteady  aerodynamic  coefficients 
along  the  span  are  presented.  Curves  for  the  wing  alone  are  drawn 
for  the  results  of  References  1 and  5 and  for  the  results  of  this 
paper.  The  curves  for  the  wing-body  combination  are  also  drawn  on 
the  same  graphs  for  comparison. 

It  is  interesting  to  note  that  excellent  agreement  between 
lift  curves  is  found  for  the  wing  alone  case  over  the  inboard  half 
of  the  span.  The  lack  of  agreement  near  the  tip  is  caused  partly 
because  of  the  small  number  of  boxes  in  this  region  and  partly  be- 
cause the  calculations  for  some  of  the  coefficients  entail  working 
with  small  differences  of  large  numbers.  The  wing -alone  moment 
curves  are  in  good  agreement  over  the  inboard  half  of  the  wing  for 
some  of  the  coefficients,  although  the  agreement  is  poor  in  other 
casea.  This  is  to  be  expected  since  the  moment  coefficients  depend 
on  the  distribution  of  pressures  rather  than  total  farces  as  in  the 
case  of  the  lift  coefficients  and  accurate  pressure  distributions 
cannot  be  expected  from  the  coarse  grid  used  in  these  calculations. 
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The  changes  in  total  lifts  and  moments  are  not  as 
extreme  as  may  be  expected  from  an  inspection  of  Figures  3 
through  10»  Large  changes  do  occur  in  the  real  components 
of  the  force  and  moment  vectors  due  to  translatory  oscillations. 
However,  since  these  real  components  are  small  compared  to  the 
imaginary  components,  the  total  changes  in  lift  and  moment 
magnitude  are  essentially  those  due  to  the  changes  of  the 
imaginary  components. 

The  overall  effect  of  the  body  on  the  wing  can  be  seen 
more  clearly  in  Figures  11  and  12.  The  ratios  of  the  magnitudes 
of  the  aerodynamic  coefficients  for  the  wing-body  combination  to 
the  wing  alone  coefficients  are  presented  in  Figure  11.  In  this 
figure,  the  following  definitions  apply. 

Lh=n 7^7 

“hi1 

It  is  seen  that  the  lift  coefficients  are  affected  most 
in  the  middle  region  of  the  wing  (the  dashed  portions  of  the 
curves  indicate  extrapolated  data),  whereas  tha  moment  coefficients 
are  altered  most  near  the  wing-body  intersection.  An  interesting 
observation  is  that  the  curves  for  the  outboard  half  of  the  semi-span 
are  almost  exactly  the  same  for  all  four  coefficients.  The 
equations  derived  in  this  paper  are  too  complex  to  determine  whether 
this  is  the  general  case  or  whether  it  is  only  a coincidence  for  this 
particular  set  of  parameters. 

Figure  12  presents  curves  showing  the  variation  of  the 
change  in  phase  angle  of  the  aerodynamic  vectors  along  the  wing 
span.  It  is  seen  tlat  for  all  coefficients,  the  body  produces  almost 
no  phase  change  over  the  outboard  two-thirds  of  the  wing.  The  lift 
coefficients  suffer  very  little  phase  change  near  the  root,  whereas 
an  extremely  large  phase  change  takes  place  near  the  root,  for  the 
moment  coefficients.  The  size  of  the  phase  shift  may  be  considerably 
in  error.  As  explained  previously,  accurate  pressure  distributions 
cannot  be  expected  when  a very  coarse  grid  is  used.  Inspection  of 
Figure  2 shows  that  the  local  chord  of  the  lattice  representing  the 
wing  changes  from  two  to  three  boxes  at  just  about  the  same  span- 
wise  station  that  the  sudden  change  in  phase  shift  occurs.  Actually 


(«7) 

(68) 

(69) 

(70) 
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the  effect  of  the  change  in  chord  is  not  as  great  as  the  above 
statement  might  imply  because  the  potential  contributed  by  each, 
box  is  weighted  according  to  the  percentage  area  of  the  box 
covered  by  the  wing.  It  should  be  noted  however  that  the  large 
changes  in  the  moment  coefficient  magnitudes  and  phase  angles 
near  the  wing-body  intersection  may  be  grossly  exaggerated. 

The  singularity  at  the  body  axis  causes  the  velocities  induced 
by  the  body  at  points  very  close  to  th  body  axis  to  be  much 
too  large,  thereby  exaggerating  the  jj-fect  of  the  body  on  the 
wing  air-loads  in  this  region. 

An  effect  which  is  not  evident  from  the  equations 
presented  in  the  previous  sections  was  noticed  during  the  calcu- 
lations. It  was  found  that  one  of  the  most  dominant  terms  involved 
in  the  calculation  of  the  body -induced  dowiwash  was  the  term 
involving  S'.  Not  only  was  this  term  larger  than  the  others  but 
it  was  noticed  that  if  this  term  were  omitted,  the  body-induced 
downwash  would  be  of  an  entirely  different  order  of  magnitude  and 
possibly  even  of  a different  sign.  Thus  if  a circular  cylinder 
were  used  as  a fuselage,  the  results  would  be  quite  different.  It 
seems  likely, therefore,  that  the  effects  of  the  body  on  the  wing 
air-loads  are  quite  sensitive  to  body  shape  and  the  position  of  the 
wing  on  the  body.  The  results  of  this  example  should  therefore 
not  be  considered  to  be  typical.  However,  they  do  show  that  wing- 
body  interference  effects  can  be  significant. 
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SECTION  V 


SUMMARY  AND  CONCLUSIONS 


An  approximate  method  has  been  developed  for  taking  into 
account  the  effect  of  a fuselage  on  the  unsteady  air-loads  an  an 
elastic  wing  with  supersonic  edges.  A numerical  procedure  for 
obtaining  solutions  is  presented,  based  on  the  technique  employed 
in  Reference  1 wherein  the  wing  alone  case  is  considered. 

The  results  of  sample  calculations  show  that  a body  with 
maximum  diameter  equal  to  twenty-five  percent  of  the  maximum  wing 
span  appreciably  affects  the  magnitude  of  the  unsteady  ae.'odynamic 
forces  and  moments,  whereas  the  phase  angles  of  the  aerodynamic 
force  and  moment  vectors  are  only  slightly  affected. 

It  became  apparent,  during  the  course  of  the  sample 
calculations,  that  downwash  induced  on  the  wing  by  the  body  is 
very  sensitive  to  the  shape  of  the  body.  Therefore  the  results 
obtained  in  the  sample  calculations  should  not  necessarily  be 
considered  as  representative  of  the  general  effect  of  a body  on 
the  unsteady  air -loads  on  a wing. 

Although  the  example  calculated  he rain  was  concerned  with 
the  effect  of  a body  on  the  airloads  of  a wing  with  supersonic 
edges,  the  method  is  also  applicable  to  a wing  with  subsonic  edges, 
since  the  essence  of  the  method  is  merely  the  modification  of  the 
do vn wash  in  the  area  perturbed  by  the  wing. 
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FIGURE  2.  CONFIGURATION  FOR  ILLUSTRATIVE  EXAMPLE 
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